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Controlling Fusion of Majorana Fermions in one-dimensional systems by Zeeman Field 
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We propose to realize Majorana fermions (MFs) on an edge of a two-dimensional topological 
insulator in the proximity with s-wave superconductors and in the presence of transverse exchange 
field h. ft is shown that there appear a pair of MFs localized at two junctions and that a reverse in 
direction of h can lead to permutation of two MFs. With decreasing h, the MF states can either be 
fused or form one Dirac fermion on the 7r-junctions, exhibiting a topological phase transition. This 
characteristic can be used to detect physical states of MFs when they are transformed into Dirac 
fermions localized on the 7r-junction. A condition of decoupling two MFs is also given. 
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In recent several years, how to realize, manipulate, and 
detect Majorana fermions (MFs) is one of the most ac- 
tive topics of research in condensed matter physics. |l|- 
5] . The nonabelian character of the MF makes it to be 
a promising candidate for topological quantum compu- 
tation [fj-y. There are many systems that manifest the 
MF, such as half-quantum vortices of p-wave supercon- 
ductors [fj, y, [lfj , the hexagonal spin lattice model [7| , 
the one-dimensional (ID) p-wave lattice [ll|, the topo- 
logical surface state with proximity to an s-wave super- 
conductor [l|, [FJ], ultracold atom systems [13J, LL4[ and 
so on. Since Majorana bound states are superpositions 
of electrons and holes in the middle of superconducting 
gap, they are neutral zero-energy excitations, and the 
particle-hole symmetry causes the antiparticle of an MF 
to be itself in the field-theory framework [l5( ■ Many pro- 
posals have been suggested to explore novel properties 
of MFs, such as the electrically detected Majorana inter- 
ferometry [3fl , the Andreev reflection induced by MFs [4| , 
the charge transport with Majorana edge mod es J 161] . and 
the teleportation by Majorana bound state [17|. For a 
superconducting system only Cooper pairs can be cre- 
ated and annihilated, and so nonabelian statistics of the 
MFs can only be formulated in subspaces of same fermion 
parity [6(. Also, manipulating MFs in ID systems can 
be achieved by using assistant quantum systems such 
as Coulomb blockaded quantum dots [18[ and semicon- 
ducting wire networks composed of trijunctions [19j. Al- 
though it has been reported recently that signature of 
experiment supports the existence of MFs [20], how to 
detect MFs still remains as an open question. 

In this paper, we propose to realize the MFs on the 
edge of a 2D topological insulator. Superconducting or- 
der parameters with amplitude A are introduced to the 
edge states by proximity effect of s-wave superconduct- 
ing junctions, and a local transverse Zeeman field h is 
also applied there. It is found that existence of a pair of 
MFs depends explicitly upon the relative magnitude of 
h and Ao |2l| . As | h |> Ao, two MFs in different spin 
components emerge at the two junctions, respectively, 



and an inverse Zeeman field will lead to permutation of 
MFs. At | h |= Ao, there will be a topological transi- 
tion. For | h |< Ao, the MFs can either be fused or form 
Dirac fermions localized on the junctions, depending on 
the phase differences of the junctions. When the phase 
difference is unequal to (2A^ + l)ir, the wavefunctions of 
MFs are extended into the bulk and fused; otherwise, 
one additional MF is created at the junction and com- 
bined with the original one to form one Dirac fermion. 
Therefore, when the phase difference of only one junction 
is equal to (27V + l)7r, the MF on the junction is effec- 
tively driven to the other junction and forms one Dirac 
fermion there. This character can be used to detect phys- 
ical states of two MFs. The coupling between two MFs 
will vanish under some conditions. In the present pro- 
posal, all the processes allow us to realize, manipulate, 
and detect the MFs readily. 
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FIG. 1: (Color online) A half-infinite 2D topological insu- 
lator (yellow region) and its ID edge (red line) along the x 
direction. Three blocks of superconductors (green blocks) are 
fabricated on its top to form two superconducting junctions at 
x — and x = L. A local Zeeman field along the y direction 
is applied to the middle (brown) region of < x < L. 



It has been shown theoretically and experimentally 
that the quantum spin Hall effect can be realized in 
HgTe/(Hg,Cd)Te quantum wells and the gapless edge 



state is protected by time reversal symmetry 22|, [23 1. 
In Fig. [TJ we consider a half-infinite 2D topological in- 
sulator, and its edge is in the proximity to three blocks 
of s-wave superconductors with different phases. The 



Hamiltonian of this system may be written as [3[ 



H 



dx{ipl l (x)(pa 1 + h{x)a 2 - fi) a pipp{x) 



free of the Zeeman field, the quasiparticle spectrum for 
plane waves is simply given by E = \/k 2 + Aq. For 
the bound state of E = 0, Hamiltonian © can be 
solved to yield (u t(i) ,u m) ) T = e ±xAo (l, ±ie~^ L ) T for 



-Aipl(x)ipUx) + A*ipi(x)i/j^(x)}. (1) x < 0. The wavefunction proportional to e 



-xA 



is m- 



Here the first term is the Hamiltonian of the ID edge 
state with a uniform Zeeman field ft, along the the y direc- 
tion for < x < L, in which p = —id x is the momentum 
operator, <7j are the spin Pauli matrices, /i is the Fermi 
energy. Also, % = Is and Fermi velocity Vp = 1 have 
been taken. Hamiltonian (fTJ) is not invariant under time 
reversal because of the presence of Zeeman field ft. The 
superconducting order parameter is given by 



A(x) = 



A e l(fiL 


x < 0, 


A 


< x < L 


A e lipR 


x > L, 



(2) 



as shown in Fig. [TJ with A > and cp's as the 
phase of each superconducting region. Obviously, 
the charge conjugation symmetry is preserved. The 
quasiparticle operator in the Nambu representation 
\^{x)} — [>pf(x),^pi(x),'>pj_(x),—')pt(x)] T is defined as 

7 = J dx{u^ilj^(x)+u* i ip ir (x) + v^i(;^(x)-v^^(x)}. When 
the quasiparticle has relation 7 = 7T, it is a neutral 
MF. The fact that quasiparticle annihilates itself leads 
to u-j- = —vl and uj_ — vT. By calculating the equa- 
tion of motion given by Ej = [7,-ff], we recover the 
form of Bogoliubov-de Gennes (BdG) Hamiltonian in 
(ut,ui,v l ,v<t) T as 

W-BdG = pcriT 3 + h(x)a 2 - M T 3 + ReA-ri - ImAr 2 . (3) 

Here r is the Pauli matrices in the Nambu representation. 
It has been pointed out that the Dirac field with the s- 
wave superconducting order parameter is equivalent to 
that of a p + ip superconductor that has zero-mode MFs 
[2J-|26[. In Eq. <j3j> , the existence of zero modes relies 
only on vanishing of the determinant for Hamiltonian 
([3|). As a result, the wave vector is readily solved as 
k = ± \/ [i 2 — ft 2 ± i Ao in the middle region. The imagi- 
nary wave vectors indicate that there are localized states 
that may give rise to the MFs. It should be noted that 
the first term of Hamiltonian (TTJ) in the absence of Ao 
breaks the time reverse symmetry and yields two Fermi 
zero modes localized at x — and x — L, which decay 
as e~\ h \ x . If we choose that superconducting order pa- 
rameters are introduced at \i = 0, only two Fermi zero 
modes have contribution to superconducting condensa- 
tion. When /! = 0, the continuous spectrum of Eq. ([3]) in 
the middle region is given by E = ^/k 2 + (h ± Ao) 2 with 
two energy gaps A± =| h ± A |. Obviously, the topo- 
logical phase transition happens when the gap is closed 
at h = ±Ao- 

We first focus on the left junction at x — in the 
case of ft- > Aq. Since the left domain in Fig. Q] is 



valid because it diverges as 1 — > — 00, so that we have 
(w-|-(_L),u-j-(4_)) T = e xA °(l,ie~ tVL ) T . For the zero-energy 
mode of x > 0, one finds that solutions are decoupled 
into the spin- up and spin-down components in Eq.©, 
yielding (u t ,v t ) T = Ae xX +(l,i) T + Be xX - (1, -i) T and 
(u i ,v i ) T = Ce~ xX +(l,-i) T + De~ xX -(l,i) T with A± = 
(h ± Aq)/ hyp. Since the solution in spin-up component 
diverges as e xX± (A± > as h > A ) for x —> 00, we have 
A = B = 0, and coefficients C and D arc determined by 
the boundary condition of wavefunction continuousness 
at x = 0. As a result, after gauging away the phase fac- 
tor, the wavefunction of E = in the left junction is 
obtained as 



IV'l) 



for x < and 



„iAo 



(0,e 



;(2£- 






,0) T (4) 



\^l) = e-* A +sin^(0,e^, e -7\0) T 

+e- xA -cos^(0,e-^,e^,0) T (5) 

for x > 0. It can be readily shown that the quasiparticle 
operator defined above as 7l = (V>l| v I'} satisfies relation 
-fL = 7lj an d so there is an MF in spin-down component 
localized at x = 0. 

The same approach can be applied to the right junction 
at x — L. It is found that there exists an MF in spin-up 
component located at X = L, whose wavefunctions are 
obtained as 



Mr) 



= e(*- L > A +sin^(e^,0,0,- e -^) T 

+e ( "- L)A -cos^(e-^,0,0,-e^) T (6) 



for x < L, and 

for x > L. Obviously, 7^ = (ipp^) satisfies jp = 7^. 
Away from each interface, the MF has two decay lengths: 
1/A + and 1/A_, which are closely related to two energy 
gaps A± =| h ± A |. Which one dominates the decay 
rate of the MF depends upon the phase difference of the 
junction. For example, if <pl — 2Nir with N an arbitrary 
integer, the left MF decays as e~ xX - for x > 0, while at 
(fL = (2N + 1)tt, it decays as e~~ xX+ for x > 0. 

In Fig. [2l we plot zero-energy quasiparticle probability 
I u \ 2 for ft, > A (a) and h < — A (b) as ipL = <pit = 0. 
As shown in Fig. 2(a), there are two MFs localized at 
x = and x = L for ft, > Ao. As the magnitude 
of the Zeeman field is turned down, the MFs become 



more and more extended. And when \h\ < Ao, the two 
MFs are annihilated on the junctions and fused. Such 
a fusion of MFs arises from the sign reverse of A_ at 
h = A (A + at h = — A ) due to the closing of en- 
ergy gap A_ (A + ). For | h |< A , since A + > and 
A_ < 0, the second terms in Eqs. (|5|) and ^ would di- 
verge as x — > oo, and so the solution for the MFs would 
be an unphysical result. Therefore, there is no MF for 
| h |< A , and there appears a topological transition at 
| h |= Ao from the MF phase to the trivial one with- 
out MF. If the direction of the Zeeman field (h > Ao) 
is reversed, we have h < — Ao and so A± < 0, with 
the result being shown in Fig. 2(b). In this case, using 
the same procedure of calculation, the wavefunctions on 
both sides of each interface at x — or x — L can be 
obtained as follow. The MF located at x — is now 
obtained in the spin-up component, and its wavefunc- 
tion is given by \i[> L ) — e xX - sin^(— e~^ l ,0, 0, ei l ) T + 



e xX + cos ^r(e4 4 ; 0, 0, — e~* 4 ) for x > 0. Compared 
with Eq. ([5]), one finds that the direction reverse 
of the exchange field leads to that the MF in the 
spin-down component is replaced by an MF in the 
spin-up component, accompanied with an exchange 
of exp(— x\±) and exp(xAzp) due to the sign reverse 
of A±. Similarly, the wavefunction of the MF lo- 
cated at x — L is obtained in the spin-down compo- 
nent as \ip' R ) = -e- ( - x - L ^-8in*f(0,e-% i ,e : z i ,0) T + 



-(x-L)X + 



COS • 



-(O.e- 



, 0) for x < L. It can also 



be obtained by performing a symmetry operation to re- 
cover the solution of a reversed exchange field. Hamil- 
tonian JTJ satisfies UH(h)W = H(-h) with U ^= KO 
where 1Z is a spin-rotation of n around x and O is the 
central inversion with x — > L — x. It is obvious that 
Wl,r) ~ u \^l,r)- 
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FIG. 2: (Color online) Space distribution of probability den- 
sity \ u \ 2 of the MFs for different parities (red and blue lines) 
with ipL = <fiR — 0- The exchange field is taken to be h — 2Ao 
(solid line) and h = l.lAo (dashed line) in (a), and h = — 2Ao 
(solid line) and h = — l.lAo (dashed line) in (b). 

The above discussion is suitable to the case of (fiL,R ^ 
(2N + 1)tt. For <p LtR = (27V + 1)tt, the second terms 
of Eqs. (0 and ([6]) vanish, and the calculated results 
for zero-energy quasiparticle distributions are plotted in 
Fig. El It is found that there exist still two MFs on the 
junctions for h > Ao. The essential difference is that 
the MF peaks have merely a very slight extension with 
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FIG. 3: (Color online) Space distribution of | u | 2 for different 
parities (red and blue lines) with ipL — <{>r = fl". The exchange 
field is taken to be h = 2Ao (solid line) and h — l.lAo (dashed 
line) in (a), h = 0.9A in (b), h = -0.9A in (c), and h = 
— 2Ao (solid line) and h — — l.lAo (dashed line) in (d). 



decreasing h, as shown in Fig. (3]Ja) and (b). At the same 
time, another pair of MFs are generated just after clos- 
ing energy gap A_ = \h — Aq| so as to form two Dirac 
fermions at x = and x = L. Such a novel behavior can 
be understood by the following argument. Taking the 
left junction for example again, the MF state in Eq. (7) 
decays as e~ xX + for x > 0. For \h\ < A , regardless of the 
closing of gap A_, A + = h + Ao makes the original MF 
survive, for this MF is protected only by gap A + . More 
interestingly, there appears an additional MF in the spin- 
up component at x = 0, whose wavefunction is propor- 
tional to e xX ~ with A_ < 0. This MF must be combined 
with the original one, forming a Dirac fermion. Since a 
change of the phase for a superconductor will not close 
energy gap of the bulk, the present phase with the Dirac 
fermion located at the junction for ifiL.R = (2JV + l)n 
is topologically equivalent to that without Dirac fermion 
there for <pl,r ¥" (2AT + 1)-7T. We wish to point out here 
that the MF in the spin-up component at x — exists 
only for Ao > h, and annihilates for h > Ao due to the 
sign reverse of A_ at h — Ao. The same argument is 
also applied to the right junction. As a result, a pair of 
Dirac fermions are formed in Fig. [3] (b). For h < 0, the 
particle distribution is shown in Figs. |3][c) and (d). With 
increasing the magnitude of h < 0, the topological phase 
transition occurs once again at h = — Ao. 

For ip L = (2iV+l)7r and ip R ^ (27V+1)7T, the situation 
is also interesting, and the space distribution of \u\ 2 with 
different Zeeman fields is plotted in Fig. |4j For h > Ao 
in Fig. Ua), the right MF extends rapidly into the bulk 
with decreasing h, whereas the left MF remains almost 
unchanged. As h is less than Ao, the MF on the right 
junction is annihilated and at the same time another MF 
is created on the left junction, as shown in Fig. [31(b). 
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FIG. 4: (Color online) Space distribution of | u | 2 for different 
parities (red and blue lines) with ipL — n and ipn = n/3. The 
exchange field is set h = 2Ao (solid line) and h — 1.1 Ao 
(dashed line) in (a), h = 0.9A in (b), h = -0.9A in (c), 
and h — — 2Ao (solid line) and h — — l.lAo (dashed line) in 
(d). 



This evolution is equivalent to the process that the MF 
on the right junction is driven to the left junction and 
two MFs there are combined to form one Dirac fermion. 
As h is reversed, the two MFs that combine into one 
Dirac fermion located at x — exchanges their magni- 
tudes [see Fig. BJc)]; and as h < — Ao, an MF moves 
back to the right junction and the system reenters the 
topological phase [see Fig. HJd)]. The underlying physics 
has been discussed above, and will not be repeated here. 
The evolution from MFs to Dirac fermions in Fig. 0]can 
be used to detect the physical state of MFs. We can drive 
them to form Dirac fermions for detection and initialize 
MFs between subspaces of different parities. 

The MFs on the right and left junctions are coupled 
with each other, i.e., {^L^BdG^RI 7= 0- F° r h > A , 
IV'l) and \ipg) have been given by Eqs. dU and ([5]), and 
© and (J7J, respectively. It can be shown that such a 
coupling depends to a great degree upon phase differences 
of the two junctions, and it will vanish if the following 
condition is satisfied, 



(PL-<P ,<P-<PR 
cot cot 



-2L/i 



(8) 



where £ = Hvf/^-o is the superconducting coherent 
length. In this case, we have a pair of zero-energy MFs 
decoupled exactly, such as those in Figs. (3) and (4) 
where ipj, = tt and/or ifR = n and L ^> £. If condi- 
tion (8) is not satisfied, the coupling will make the MFs 
have a small departure from zero energy, proportional to 

In summary, we have shown that the edge state of a 2D 
topological insulator in the proximity with s-wave super- 
conductors and under a vertical Zeeman field may accom- 



modate MFs. The MFs can be manipulated by tuning 
the phase differences of the junction to either be fused 
or form one Dirac fermion on the 7r-junction. As the ex- 
change field becomes less than the superconducting order 
parameter, one MF can be driven from a non-7r-junction 
to a 7r-junction, forming one Dirac fermion on the n- 
junction. There exhibits a topological phase transition 
at I h |= Ao. For | h \< Ao, there are three topologically 
equivalent phases with different phase differences: one 
and two Dirac fermions localized on the 7r-j unctions, and 
none. We have also illustrated a condition that decou- 
ples completely the two MFs at two junctions regardless 
of the magnitude of the exchange field. 
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